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TIME-REVERSAL AND ELLIPTIC BOUNDARY 
VALUE PROBLEMS 

By Zhen-Qing Chen 1 and Tusheng Zhang 

University of Washington and University of Manchester 

In this paper, we prove that there exists a unique, bounded con- 
tinuous weak solution to the Dirichlet boundary value problem for 
a general class of second-order elliptic operators with singular coeffi- 
cients, which does not necessarily have the maximum principle. Our 
method is probabilistic. The time reversal of symmetric Markov pro- 
cesses and the theory of Dirichlet forms play a crucial role in our 
approach. 

1. Introduction. The pioneering work by Kakutani [18] in 1944 on rep- 
resenting the solution to the classical Dirichlet boundary value problem 

Au = 0, in D, 
u = f, on dD, 

using Brownian motion started a new era in the very fruitful interplay be- 
tween probability theory and analysis. Here D is a bounded connected open 
subset of M. n . Since then, in place of Laplacian A, there are two classes of 
second-order elliptic differential operators that have been studied in connec- 
tion with probabilistic approach. One is the nondivergence form operator 

1 n n d 

(L1) £= 2^^ W ^ + ^' W &- + ^ 

1,3=1 ■> i=l 

The other is the divergence form operator 

1,3=1 i=l 
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where A{x) = (a,ij(x)) is an n x n symmetric bounded positive definite ma- 
trix. For nondivergence form operator C in (1.1), one can run stochastic 
differential equation 

dX t = a(X t )dB t + b(X t )dt, 

where a is a symmetric nxn matrix such that a 2 = A, b = {b\ , . . . , b n ) and B 
is a Brownian motion on M 71 . The infinitesimal generator of X is Cq = C — q. 
Under some suitable conditions, the solution for Cu = in D with u = f on 
dD can be solved by 



u(x) = E 



e^J T \{X s )ds^f{X TD ) 



for x € D; 



see [12]. Here td = inf{t > 0:X t ^ D} is the first exit time from D by X. 
When C is the divergence form operator of the form (1.2), one has to run 
symmetric diffusion associated with ^V(^4V). Observe that X is in general 
not a semimartingale when A is just measurable. Nevertheless one can still 
use the symmetric diffusion X to solve the Dirichlet boundary value problem 

Cu = in D with u = f on dD 

through a combination of Girsanov and Feynman-Kac transforms (see [8]). 
The Dirichlet form theory plays the role of Ito's calculus in the divergence 
form operator case. 

In this paper, we study the Dirichlet boundary value problems for second- 
order elliptic operators of the following form: 



(1.3) 



C = -V • (AV) + b • V - div(b-) + q 

1,3=1 i=l 



in a bounded domain L> C R n . Here A = (%■) :M n — > M n x ]R n is a Borel 
measurable, symmetric matrix-valued function that is uniformly elliptic and 
bounded, that is, there is a constant A > 1 such that 

(1.4) A _1 / nX n < A(x) < Mnxn for every x G M n ; 

b = (6i, . . . , 6 n ) and b = (6i, . . . , b n ) are Borel measurable M n - valued func- 
tions on M n and q is a Borel measurable function on IR n such that 

(1.5) l D (|b| 2 + |b| 2 + |g|)GK n . 

Here K n denotes the space of Kato class of measures on R n : when n > 3, a 
signed measure (i £ K n if and only if 



lim sup / \x — y\ 2 n \fi\(dy)=0, 

r ^°xm n JB(x,r) 
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where \p,\ denotes the total variational measure of \i. Kato class K n can 
also be defined for n = 1,2; see [9] for details. A function q is said to be in 
K n if fi(dx) := q{x)dx is in K n . Clearly L°°(R n ) C K n , and it is easy to 
see by using Holder's inequality that L p (M. n ) C K n for p > n/2. By taking 
b = 6 = = goff-D, we may and do assume in the sequel that (1.5) holds 
without the restriction of D by Id- 

Note that "div(b-)" in (1.3) is just a formal writing because the divergence 
really does not exist for the merely measurable vector field b. It should be 
interpreted in the distributional sense. It is exactly due to the nondifferentia- 
bility of b, all the previous known methods in solving the elliptic boundary 
value problems such as those in [8] and [12] ceased to work. The lower-order 
term div(b-) cannot be handled by Girsanov transform or Feynman-Kac 
transform. We will show in this paper that this term in fact can be tackled 
by the time-reversal of Girsanov transform from the first exit time td from 
D by the symmetric diffusion X associated with ^V(-AV). This is the nov- 
elty of this paper. Note that time reversal of a Girsanov transform from a 
deterministic time was first studied in [20] for diffusions, and very recently 
in [4] in the context of general m-symmetric Markov processes. We point out 
that time reversal from a deterministic time in [4, 11, 20] are defined under 
the stationary measure P m . Doing time reversal from a random time td in- 
volves many delicate technical issues for an effective analysis. We are able to 
circumvent these difficulties through a certain /i-transform. See (1.12) below 
for details. 

Let (Q,T>(Q)) be the bilinear form associated with the operator £, where 



For an open subset D C R", denote by C£°(D) the space of smooth functions 
on D with compact support. The L 2 -domain T>(£) of £ is defined to be 



and we denote g by £u. Clearly, it follows from this definition that 



V{Q) = {u G L 2 {R n ) : Vn G L 2 (R n )} = W 1 ' 2 (M") 
and for u, v G W 1,2 (R n ) 




2(«.^) = o E / a ij( x )x— dx -Y\ / bi(x)—v{x) 
2^ 1 7k« dxidxj ~[J«- n dxi 



{u G W 1,2 (K n ): there is g G L 2 (W 



) so that Q(u,v) = {-g,v) L 2 {Rn) 
for every v G C^°(IR n )} 



Q(u,v) = (—£u,v) L 2(^n) for u G T>(£) and v G T>(Q). 
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It is well known that the differential operator £ enjoys the maximum 
principle if — div(b) + q < in W l in the following distributional sense: 

(1.7) V/ bi(x)p-dx+ [ q(x)<j)(x)dx<0 

JR" OXi Jr« 

for all nonnegative function <f> in C£°(M n ). 

Triidinger [26], Theorem 3.2 and Corollary 5.5, has proved the following: 

Theorem 1.1. Assume the Markovian condition (1.7) holds. For every 
f £ W > 2 (D), there exists a unique weak solution u 6 W 1,2 (D) such that 

(1.8) £u = in D with u — f £ Wq' 2 (D). 

Here Wq' 2 (D) is the completion of C£°(D) under the Sobolev norm 




in W 1,2 (D). Moreover u is locally Holder continuous in D. 

Recall that Cu = in D is understood in the following distributional sense: 

Q(u,<P)=0 for every € C7 C °° (D) . 

We stress that condition (1.7) plays a key role in Triidinger's approach be- 
cause of the critical use of the maximum principle there. 

The aim of the present paper is twofold. The first is to give a proba- 
bilistic representation for the weak solution of the Dirichlet boundary value 
problem (1.8). This is highly nontrivial because there is no longer a Markov 
process associated with the operator £ due to the appearance of the lower- 
order term div(b-), nor can that lower-order term be handled via Girsanov 
transform or Feynman-Kac transform. Our idea is to use the symmetric dif- 
fusion process X associated with the divergence form operator ^V(AV), the 
symmetric part of £, and treat £ as its lower-order perturbation via a com- 
bination of Girsanov and Feynman-Kac transforms and a time-reversal of 
Girsanov transform at the first exit time td from D by X. Based on the new 
probabilistic representation, our second aim is to establish the existence and 
uniqueness of the weak solution to problem (1.8) without the Markovian as- 
sumption (1.7). To this end, we introduce a kind of /i-transformation which 
transforms the solution of the problem (1.8) to the solution of a Dirichlet 
boundary value problem for operators which do not involve the adjoint vec- 
tor field like b. The time reversal and the theory of Dirichlet forms play an 
essential role throughout this paper. 

The remaining of the paper is organized as follows. Let X be the symmet- 
ric diffusion with infinitesimal gener ator \V{AV). It is well known (cf. [25]) 
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that X is a conservative Feller process on W 1 that has Holder continuous 
transition density function which admits a two-sided Aronson's Gaussian 
type estimate. In general, X is a not a semimartingale but it admits the 
following Fukushima's decomposition (cf. [14]): 

(1.9) X t = X + M t + N t , t>0, 

where M = (M, . . . , M n ) is a martingale additive functional (MAF) of X 
with quadratic co- variation (M\M^) t = / * a^(X s ) ds and N = (N 1 ,..., N n ) 
is a continuous additive functional (CAF) of X locally of zero quadratic 
variations. Note that, since X has a continuous density function and each 
dij is bounded on R d , by [13], Theorem 2, M and N can be refined to be 
CAFs of X in the strict sense without exceptional set and (1.9) holds under 
Pa; for every x € K n . Without loss of generality, we work on the canonical 
continuous path space C([0, oo), W 1 ) of X and, for t > 0, denote by r t the 
reverse operator of X from time t. In Section 2, we prove that, under the 
Markovian condition (1.7), the solution u of the problem (1.8) admits the 
following probabilistic representation: 

(1.10) u(x) = B x [Z TD f(X TD )} forxGD, 

where td := inf {t >0:Xt£ D} is the first exit time from D by the symmetric 
diffusion X and 

Z t = exJ [\A- 1 h)(X s )dM s + ( f (A-^X^dMs) or, 

(1.11) U ° t U ° ^ J 

- 1 jf*(b - b)^~ 1 (b - b)*(Jf«) + £° q(X s ) ds 

All the vectors in this paper are row vectors and we use b* to denote the 
transpose of a vector b. For two vectors a and (3 in M n , we use a- (3 or (a, /3) 
to denote their inner product. 

Note that Z TD is well defined under P x for quasi-every x £ M n . This is 
because 1 1— > /Q(74 _1 b)(X s ) dM s is a square integrable martingale of X in 
the strict sense having finite energy and so it follows from [3] and [23] that 
there is a continuous additive functional L of X having zero energy (which 
in general may admit an exceptional set) such that for quasi-every (q.e.) 
x e R d , P z -a.s. 

^V _1 b)PQ dM^j or t = - ^'(A- 1 b)(X s ) dM s + L t , t> 0. 

To prove the probabilistic representation (1.10), we introduce a sequence 
of approximating operators with b in the definition of C replaced by 
smooth bfc. We first show that the solution Uk of the Dirichlet boundary 
value problem for the operator has the representation (1.10) with b 
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replaced by smooth b&. To complete the proof, we then show that both 
the approximating solutions and the corresponding probabilistic expressions 
converge in an appropriate sense. The purpose of Section 3 is to show that the 
weak solution u to the Dirichlet boundary value problem (1.8) is continuous 
up to the boundary of the domain D. The crucial step is to show that for 
every L p -integrable i? n -valued function / in a ball Br with radius R > and 
p> n, there exists a function v G W ' p (Br) such that the following identity 
holds: 

( T f(X s ) dM s ) o rt = - [ f(X s ) dM s + N? 
(1.12) U ° 7 J ° 

for t < inf {s : X s <£ Br}, 

where N v is the zero-energy part of the Fukushima's decomposition for 
v(X t ) — v{Xq) (see Section 3). By Sobolev embedding theorem, if we extend 
v to W 1 by taking v = on B R , then v is continuous on W 1 . It then follows 
from [13], Theorem 1, that N v can be refined to be a CAF of X in the strict 
sense. Consequently, t \— > (/J f(X s ) dM s ) o n can be refined to be a CAF of 
X in the strict sense by using the expression on the right-hand side of (1.12). 
In particular, Z TD is then well defined under P x for every x € D. Equation 
(1.12) allows us to get rid of the reverse operator rt in the expression of 
the solution u making the analysis possible. In Section 4, we establish a 
general existence and uniqueness result for the weak solution of the Dirichlet 
boundary value problem (1.8) for the operator C without the Markovian 
assumption (1.7). Let Z be defined by (1.11). We establish as Theorem 4.4 
the following important gauge theorem under suitable condition of D, A, b, 
b and q, which is of independent interest 

if E Xo [Z TD ] < oo for some xo € D, then x \— > E X [Z TD ] is bounded 
between two positive constants. 

We then show that if ~E Xo [Z TD ] < oo for some xq £ D, then for every / € 
C(dD), the equation Cu = in D has a unique weak solution that is contin- 
uous on D such that u = f on 3D. Moreover, this solution can be expresses 
as 

u(x) = E x [Z TD f(X TD )] for every x G D. 

Our strategy is, by using (1.12), to reduce the Dirichlet boundary value 
problem Cu = in D with u = f on dD to a corresponding Dirichlet bound- 
ary value problem for an operator that does not have the lower-order term 
div(b-). In Section 5, we consider the special case of C in (1.3) where 
b = b = -AVp where p G W 1 ' 2 ^) with Vp G L p (R n ; dx) for some p > n. 
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By Sobolev embedding theorem, p is continuous on M n . In this case, the 
quadratic form (Q, W 1,2 (W 1 )) in (1.6) takes the following form: 

. 1 f , s du dv , 1 f , .d(uv) dp , 
2 *7± x Jru dxi dxj 2 Jru dxi 8Xj 

+ / u(x)v(x)q(x) dx 



for u,v £ W 1 ' 2 (W l ). For this case, we can establish a stronger result without 
additional condition on the diffusion matrix A. Let X be the symmetric 
diffusion with infinitesimal generator ^V(^4V), and recall the Fukushima's 
decomposition in the strict sense (cf. [13], Theorem 1): 

p(X t )-p(X ) = M t p + Nf, t>0, 

where M p is the MAF of X in the strict sense having finite energy and N p 
is the continuous additive functional of X in the strict sense having zero 
energy. Define 

Z t = exp (n£ + jT* q (X s ) ds^j, t> 0. 

We prove in Theorem 5.1 that if D is a bounded Lipschitz domain and if 
E XQ [Z TD ] < oo for some xq € D, then x i— >• E X [Z TD ] is bounded between two 
positive constants. Moreover, assuming that E XQ [Z TD ] < oo for some xq £ D, 
we show that for every / G C(dD), 

u(x):=E x [Z TD f(X TD )}, xeD, 

gives the unique weak solution of Cu = in D that is continuous on D with 
u = f on 3D. 

In this paper, we use ":=" as a way of definition. A statement is said to 
hold quasi-everywhere (q.e.) on some set A C M. n if there is an exceptional 
set N of zero capacity so that the statement holds on A\N. For the general 
theory of Dirichlet forms and Markov processes and their terminology, we 
refer readers to [14] and [21] . 

2. Probabilistic representation. In this section, we will give a probabilis- 
tic representation of the weak solutions of Dirichlet boundary value prob- 
lems. Consider the following regular Dirichlet form on W 1 : 



1 :J 

(2.1) 

V{£) = W 1 ' 2 (R n 



If du dv 

£( u , v ) = n/Z / a v( x )^—^—dx, 
2 Jm" axi dxj 
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It is well known that there is a symmetric conservative diffusion process 
X = {Xt,9t, rt,P x ,x £ M. n } associated with it. Since X has Holder continuous 
transition density function with respect to the Lebesgue measure on R n (cf. 
[25]), X can be modified to start from every point in M n . Without loss of 
generality, we may and do assume that X is defined on the canonical sample 
space il = C([0, oo) — » W 1 ) on which the time-shift operators {6t, t > 0} and 
time-reversal operators {rt,t > 0} are well defined: for t > 0, 

(2.2) t (u)(s)=u(t + s) fors>0 
and 

(2.3) «° e r t ^ 

Let {J-t,t > 0} be the minimal augmented filtration generated by the diffu- 
sion process X. For every u £ D{£), the following Fukushima's decomposi- 
tion holds: for q.e. i£l n , 

(2.4) u(X t )-u(X ) = M? + N?, Px-a.s., 

where M u a continuous MAF of X having finite energy and iV" is a CAF of 
X having zero energy. Note that the MAF M u and the CAF N u typically 
admit an exceptional set N of zero capacity in their definition. However 
since X has a Holder continuous transition density function, by [13], Theo- 
rem 1, both M u , N u and the above Fukushima decomposition (2.4) can be 
strengthened to admit no exceptional set [so in particular, (2.4) holds for 
every x € M. n ] if u is continuous and the energy measure for M u 



^ , du(x) du{x) 



H {u) {dx) := a id( x )-gZ. — g-— dx 
i,j=l 1 o 

is a smooth measure in the strict sense. The latter means that there is an 
increasing sequence of finely open sets {D^^k > 1} so that Ufc^i Dk =R n , 
Id^A*^) is a finite Borel measure and G\(lD k H( u )) 1S bounded for every k > 1. 
In the sequel we call an additive functional strict if it admits no exceptional 
set and call the refined decomposition of (2.4) without exceptional set a strict 
Fukushima decomposition. In fact by [13], Theorem 2, both M u and N u 
can be taken to be strict AFs of X and the strict Fukushima decomposition 
holds for every continuous function u that is locally in T such that 
is a smooth measure in the strict sense. Applying the above to coordinate 
functions fj(x) := Xj for j = 1, . . . ,n, we have 

(2.5) X t = x + M t + N u Pa.-a.s- 

for every x G R n , where M t = {M}, . . . , M") is a continuous local MAF of X 
in the strict sense with 

{M\Mi) t = f aij (X s )ds 
Jo 
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and Nt is a CAF of X locally of zero energy in the strict sense. In particular, 
there is a Brownian motion B = (B 1 , . . . , B n ), which is a martingale AF of 
X in the strict sense, such that 



M= f a{X s )dB s , t> 
Jo 



0. 



where o~(x) is the positive definite symmetric square root of the matrix A{x). 
Note that M is a MAF of X in the strict sense. 
Here is the first representation result. 

Lemma 2.1. Assume condition (1.7) holds and that b is C 1 -smooth. 
Then for every f G W 1,2 (D) n C(D), the unique weak solution u to (1.8) 
admits the following representation: for x G D 

u(x)=E x /(X T Jexpj^ T£ \A~ l h){X s )dM s 

+ ( f TD (A- 1 b)(X s )dM 9 

(2.6) \ J \ D 

-- f D (b-b)A~ 1 (h-hY(X s )ds 
2 Jo 

fTD "I " 

+ / q(X s )ds 
Jo 

Moreover, u G C{D). 

Proof. First we note that by [20], (46), we for a.e. x G D, P^-a.s. have 

( f TD (A- 1 h)(X s )dM s 
(2-7) U ° 



r {A~ i h){x s )dM s - r dw(b)(x s )ds. 

JO JO 



Note that since b G C^R"), both t h-> / * (^~ 1 b)(A" s ) dM s and t h-> f* div(b) x 
(X s ) ds are continuous AFs of X in the strict sense. Therefore (/ ro (^4 _1 b) x 
r TD can be refined using the right-hand side of (2.7) so that it 
is well defined under P x for every x G D. In particular, the right-hand side 
of (2.6) is well defined for every x G D. 

Under our assumptions, the operator C can be written as 



i-.r- 



n d 

Y^(bi(x) - bi(x))— - divb(x) + q(x). 

A — 1 2 
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Define a family of measures {Q x , % € M. 71 } on Too by 
where 

H t = exp Qf* A' 1 (b - b) (X s ) dM s - i jT (b - b) A" 1 (b - b)*(X s ) ds^j . 

(2.9) 

It is known (cf. [8]) that under measure {Q x , x G M n }, X is a diffusion process 
on ]R n having infinitesimal generator 

2 .4< 



1 " <9 / 9 



i=l 



By Theorem 5.11 in [8], when / e W X ' 2 {D) C\C(D), the unique weak solution 
u to (1.8) is continuous on D and has the following probabilistic represen- 
tation 



u 



(x) = E? 



/(X TD )exp 



div(b)(X s ) + q(X s ))ds 



xeD, 



where E^ stands for the expectation with respect to the measure Q x . Since 
{Ht/\ TD , t > 0} is a uniformly integrable martingale under P x for every x £ D, 
due to the fact that |b — b| 2 G K n (see [6], page 746), we have 



u(x) = E a 



f(X TD ) exp QT° (A" 1 (b - b)) (X s ) dM s 

- \ (b - b)A~\b - b)*{X s ) ds 

xexp^ D (-dw(b)(X s )+q(X s ))ds 
This together with (2.7) implies (2.6). □ 
Put 

J{x) 



Co exp 
0, 



1 - x 



if |x| < 1, 
if \x\ > 1, 



where Co > is a normalizing constant so that 

J(x) dx = 1. 

For any positive integer k > 1, we set 

J fc (x) := k n J(kx) for x S : 
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and 

bfe(x) = J k * b(x) := / h(y) J k (x - y) dy, 

(2-10) 

q k (x) = J k * q(x) := / q(y)J k (x - y) dy. 



Recall that we assume that b, b and q are set to be zero off D, and 
|b| 2 + |b| 2 + MeK n . 

LEMMA 2.2. (i) b k — > b in Lf oc (R n ) and q k — > q in L^ c (R n ) as k —> oo. 

(ii) For every nonnegative function <f> in C^°(M. n ) and k>l, 

Y] [ hi(x)ir-dx+ [ q k (x)4>(x) dx < 0. 
Jm.n dxi jRn 

(iii) Assume n>3. Then 

r f H{y)\ 2 + \g k {y)\ 

hm sup / : —75 dy = 0. 

k>l J{ye«. n : \x-y\<r} \X - y\ n 
x£R n 



(iv) Assume n > 3. Then 

hm sup / \^y)-^y)\ 2 + \^y)-^y)\ dy=0 , 

k^o°x€_K n JD \x-y\ n 



Proof, (i) follows easily from the fact that b G Lf oc (R n ) and q G L 1 1 0C (]R n ) 
For a nonnegative function <j> in C£°(M n ) and z G M™, put cp z (x) = 4>(x + 2). 
Then it is easy to see that for k > 1, 



Y] / 6fci(x)— / 

jRn dXi jRn 



q k (x)<f>(x)dx 



I Jk(z)(it I b i (x)^ ± dx+ ( q( 
jR n \'~ 1 JR n oxi J«n 



x)(f> z (x) dx dz 



and so (ii) follows from (1.7). Since |b| 2 G K n , (iii) is a consequence of the 
following inequality: 



|b fc | 2 (y) d < f j Rn J k (z)\b\ 2 (y-z)dz 

x&g"J\x-y\<r \x - y| n ~ 2 ~ J\x-y\<r \x - y\ n ~ 2 

k>l 



b\ 2 (y- z) 

\x-y\<r \x~y\' 



Mz)[l. . 'J V^ l'fe 
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\J\x-z-y\<r \X - Z - y\ n 1 J 

f \b\ 2 (y) . 
< sup / -j j— ^ ay, 

xGR n J\x—y\<r \X - y\ n 1 

where we used the fact J Rn Jk(z) dz = 1. The proof for q is similar. To prove 
(iv) we observe from the proof of (iii) that for any r > 0, 

f |b fc -b| 2 (^ 
sup / — ay 

<2sup / | i»- 2 d y+-^2 f \bk-h\ 2 (y)dy. 

xGM. n J \x-y\<r \x - y\ n z r n A Jr« 

Using (i), this implies that 

f \h k -b\ 2 (y) , f \b\ 2 (y) 
lim sup / — 9 ay < 2 sup / -j j — ~ ay 

ISM" JK" F-2/1" i6K n ^|ic-i/|<r F-2/r 

for any r > 0. Letting r — > we get (iv) for b& — b. The proof for \q k — q\ 
goes in a similar way. □ 

For integer k > 1, define 

1,7 = 1 1 x J 

(2.11) 

+ (b(x) - b fc (x)) • V - divb fc (x) + q k {x), 

where h k ,q k are the functions defined by (2.10). Denote by Qfc(v) the 
quadratic form associated with Ck- _ 

Now we can drop the assumption of b being C 1 from Lemma 2.1 by using 
the smooth approximation h k for b. 

Theorem 2.3. Suppose that condition (1.7) holds and f € W 1,2 (D) n 
C{D). Then the unique weak solution u of (1.8) has the following probabilis- 
tic representation: for q.e. x £ D, 



u(x) = E a 



+ 

(2.12) 



f{X TD ) exp^V^b)^) dM s 

J^ D (A- 1 b)(X s )dM s ^ or TD 

f TD (h - b)^~ 1 (b - b)*{X 8 ) ds 
Jo 

+ / q(X s )ds 
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Proof. For simplicity, we assume n > 3. [The case of n = 1 and n = 2 
can be handled similarly with a corresponding version for Lemma 2.2 (iii) and 
(iv).] Recall the differential operator C k defined by (2.11). Let u k denote the 
unique weak solution of the following Dirichlet boundary value problem: 



£ k u k = 0mD with u k - f e Wl' 2 {D). 
Since h k is smooth, it follows from Lemmas 2.2(h) and 2.1 that 

u k (x) = B x 



f(X TD ) exp ( jT" (A^b)(X s ) dM s 

(£ D (A- l b k )(X s )dM s 

1 f TD (h-h k )A~ 1 (b-h k y(x s )d s 

2 Jo 



+ 



+ 



q k {X s )ds 



Let v denote the right-hand side of (2.12). We will show that lim^oo u k (x) 
v(x). To this end, put 



Z k := exp^ TD (,4- 1 b)(X s ) dM s + (^fj {A- l h k )(X s ) dM t 

- \ (b - h k )A-\h - hkT(Xs) ds + £° q k (X s ) ds^j 



and 



Z := exp (J (A- l h)(X s ) dM s + (J (A^h^X,) dM s \ o r T 

- (b-b)A- 1 (h-h)*(X s )ds+ q(X s )ds 

2 J o JO 

We first prove that Z k — > Z in probability as A; — ► oo. It is clear that 

/ (b-b fc )A- 1 (b-b fc )*(* s )ds+ / %(X s )ds 
j o Jo 

converges in probability under P x for every x € D to 

/ (b-b)A- 1 (b-b)*(X s )ds+ / g(X s )ds. 

JO JO 



Thus, it is sufficient to show that 



(2.13) (A~ l h k )(X s )dM t 



o r 



{A~ 1 h)(X s )dM 8 ) or. 



TD 
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in probability under P x as k — > oo for q.e. x £ D. Define 



Mt 



:= f{A~ 1 h k ){X s )dM s , t>0 
Jo 



and 



M t := f t (A- 1 h){X s )dM s , t>0, 
Jo 



which are MAFs of X in the strict sense of finite energy (recall that we as- 
sumed b = b = off D). It follows from [3] and [23] that there are continuous 
processes iV t fc and N t of zero energy such that 

or t = —M-t + N t and M t or t = -M t + N t . 

Moreover, since M k — > M as k — > oo with respect to the energy norm in the 
martingale space, for every subsequence {n^}, there is a sub-subsequence 

{ n kj} so that N t 3 converges to Nt uniformly on compact intervals P x -a.s. 
for q.e. (cf. [14]). Thus, for any T > 0, on {uj:td(lv) < T}, it holds 

that 

^ D (A- 1 b k ){X s )dM s ^j or TD - (^£ D (A- 1 b)(X s )dM s \ or Tfl 

< sup |M t fe o r t - M t o r t \ 

0<t<T 

< sup \M^-M t \+ sup \Nj°-N t \. 

0<t<T 0<u<T 

This proves (2.13). So, to show that u n — > v, it suffices to prove that the 
family {Z n ,n > 1} is uniformly integrable under for q.e. x € D. In view 
of Lemmas 2.1, (2.7) and 2.2(h), we have 

Z k <L k := expQ^GA-^b -b k ))(X s ) dM s 

- \ f\h - h k )A-\h - h k )*(X s ) ds^j . 

Define L in the same way as L k but with b in place of b&. Then 

f\L k -L\dP x =f (L k — L) d~P x + / (L-L k )dV x 

Jn J{L k >L} J{L k <L} 

= [ L k dP x - I LdP x + 2 f (L-L k )dP x 
Jn Jn J{L k <L} 

= 2 [ (L-L k )dP x ^0 as/c^oo. 

J{L k <L} 
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Here we have used the fact that E X [L&] = 1 = E X [L], which is a consequence 
of the Kato class assumption on |b| 2 + |b| 2 (cf. [2]). This particularly implies 
that {L k ,k > 1} is uniformly integrable under P x for every x G D, so is 
{Z k ,k>l}. 

To show that the weak solution u of (1.8) is equal to v, by the uniqueness, 
it suffices to show that v is a weak solution to (1.8). By Theorem 3.2 (and 
its proof) of Triidinger [26], there is a constant C > 0, independent of k, 
such that 

||«fc||i,2 < C||/Ili,2 for every k>l. 
By taking a subsequence if necessary, we may assume that u n converges 
weakly to some v\ in W l ' 2 (D) and that its Cesaro mean {k~ l J2j=i Uj,k> 1} 
converges to some V2 in (W 1,2 (D), \\ ■ Hi^)- Clearly v\ = V2 = v. Moreover, 
since u k — f G Wq' 2 (D), we have v — f G Wq ,2 (D). 

As Q k (u k ,<f>) = for any <f> G C£°(D), if we can show that for any eft G 
C%°(D), Q(v,4>) = lim^oo Qk{uk,<f>), then v G W 1,2 (D) is a weak solution 
for Cv = in D with v-fe Wq' 2 (D). 



Observe that 

= — ^ f n^JrA — —rlnr — \^ ( bi{x) — - (f) dx 



Qfc(u fc ,0) = -V / a ij (x)^ — dx -J2 [ 
2 Jmn dxi dxn f-f Jm™ 

~Y1 ( b ki (x)^-u k (x)dx - [ 

J]gn OXi jRn 



q k {x)u k {x)<j)dx. 



Obviously 

1 f < \ du k d< t> j 1 f i \ dv d 4> , 

I — > / (j..-.-ItI dnr ■ = — > / aij(X)— — 7: dx 

OXi OXj 



v 1 f , ,du k d(f> 1 ^ f 

hm - > / aij(x)— — — — dx = - > / 
2 frL 7 R n 13 v > dxi dxj 2 ^ Jm 



and 



Vm£[ Hx) d -^ m dx = ±( U*) 9V{X) 



dxi 

On the other hand, for (j) G C£°(D), by Lemma 2.2 we have 

d(j){x) ST 1 f T / \d ( t ) ( x ) 



4>(x) dx. 



and 



lim V / b ki (x)— u k (x)dx = V / &i(x)^^v(;c) cfa; 



lim / qk{x)uk{x)4>{x) dx = I q(x)v(x)(p(x) dx. 

k— »00 JTRn 



This proves that 

and so Cv = in D in the distributional sense. This proves the theorem. □ 



Q(v, 4>) = lim Q k (u k , (j)) = for every </> G C™(D) 

k— >oo 
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3. Continuity at the boundary. In this section, we study the regularity of 
the solution of the Dirichlet boundary value problems (1.8) at the boundary 
of the domain. First, we prepare two useful lemmas. The next result is due 
to Meyers [22], Theorem 1. 

Lemma 3.1. For every xq G W 1 , R > and p > n, there is a constant 
e G (0,1), depending only on n,R and p, such that if 

(3.1) (1 - e)I nxn < A(x) < I nXn for a.e. x G B R := B(x ,R), 

then ^V(j4Vu) = divf in Br has a unique weak solution in Wq' p (Br) for 
every f = (/i, . . ., f n ) G L 2 (Br; dx). Moreover, there is a constant c> in- 
dependent of f such that 

\\^ u \\LP{B R ;dx) < c||f \\LP(B R ;dx)- 

Observe that since u G Wq' p (D) C M. n with p > n, by the classical Sobolev 
embedding theorem (see, e.g., [16], Theorem 7.10) u G C(D) if we take u = 
on D c . Recall that D is a bounded domain in M n . Select xq G W 1 and R > so 
that B(xq,R) d D. For simplicity, denote B{xq,R) by Br- Let X Br denote 
the symmetric diffusion in Br associated with the infinitesimal generator 
^V(^4V). The time-reversal operator for X Br will still be denoted as rt- 

For a MAF M of X Br of finite energy, let 

T{M) t = -\{M t + M t or t ) ioit<T BR . 

According to the representation theorem of martingale additive functionals 
in [14] , there is a measurable function F : Br — > W 1 such that 

M t = f F{X s )dM s for t<r BR . 
Jo 

We have the following result. 

Lemma 3.2. Suppose that p> n and the diffusion matrix A satisfies the 
condition (3.1). Then for every LP -integrable W 1 -valued function F:Br — > 
W 1 , there exists a function v G W ' p (Br) C W ' (Br) with \\Vv\\ LP (B R ;dx) — 
c\\ F \\Lp(B R -,dx) such that 

(3.2) T(M) t = T(M v ) t = Nl' fort<r BR . 

Moreover, if we extend v to W 1 by taking v = on B^, then v G C(lR n ). 

Proof. By [10], Corollary 3.2, the following orthogonal decomposition 
holds with respect to the inner product induced by the energy norm. 



(3.3) 



M t = M? + K t for t < t b 
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where v is an element in W ' (B R ) , K is a MAF of X B R satisfying T(K) = 
and ti(K,M<t>)(BR) = for all <fi G D(£). Here i s the signed Revuz 

measure of CAF (K,M^) of X of finite variation. Hence, 

r(M) t = r(M*) t = N? for t < t Br . 
By the representation of MAFs, we have 

M v t = [ t Vv(X s )dM s and K t = f G(X s )dM s fort<r Bfl , 
Jo Jo 

for some measurable vector field G = (G\(x), . . . , G n (x)) : .Br — ► R n . Since 
^(k,m<?>>(- b r) = for an G Wo 1 ' 2 ^^)' we nave 




for all 4> G Cl(B R ). This says that div(^G) = in Sr. Note that 

F(a?) = V«(s) +G(x). 

Multiplying both sides of the above equation by the matrix function A(x), 
we see that v G W Q ' (Br) solves the equation 

div(AVv) = div(AF) in B R . 

By Lemma 3.1, v G Wq' p (B r ) with 

ll w l|LP(Bji;di) ^ c ll^llLP(B J j;da:) < c l l|-Fl|.W(.Bii;dx)- 

If we take v = on 1?^, then v G W ,p (R n ) and so by the Sobolev embedding 
theorem [16], Theorem 7.10, v G C(W l ). This completes the proof of the lem- 
ma. □ 

Henceforth, we select and fix a ball B R D D. 

Theorem 3.3. Assume the Markovian condition (1.7) is satisfied, p> n 
and that the diffusion matrix A satisfies the condition (3.1) of Lemma 3. 1 . 
Assume that |b| G L p (D;dx) , and that l£>(|b| 2 + \q\) G K n . Let u be the 
unique weak solution of the Dirichlet boundary value problem (1.8). Then 
for y G dD that is regular for (^A,D), we have 

lim u(x) = f(y). 
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Proof. It is enough to prove the theorem for nonnegative function /. 
Let {$7, J 7 , X t , Q x , x E W 1 } denote the diffusion process defined as in (2.8). 
As in the proof of Theorem 2.3, put 

— r tAr B R , ~ 
M t = (A^h)(X s )dM s . 
Jo 

By Lemma 3.2, there exits a bounded, function v € W ' p (Bji) C Wq ,2 (Br) 
such that 

M t ar t = -M t + N% tert<T BR . 

Note that M is a MAF of X in the strict sense and N v is a CAF of X in the 
strict sense of zero energy in view of [13], Theorem 1, since v G Wq 1,p (Bji) 
and so it is continuous on W 1 if we extend v to take value on B^. Therefore 

Mt o rt can be refined to be a CAF of X in the strict sense. It follows that 
u can be expressed as 

u(x) = B<?[f(X TD )exp(A TD )], 

where 

A t = Nf+ fq(X s )ds. 
Jo 

Under condition (1.7), the CAF A t is negative and decreasing in t G [0, tb r ). 
Hence for x € D, 

\u(x) - E$[f(X TD )}\ < H/IUE^tlexp^) - 1|] < ||/||ocEQ[l - exp(A r J]. 
We know from [8] that for y G dD that is regular for (^A,D), 

lim E$[f(X TD )]=f(y) 

x^y,x£L> 

and lim x ->y tX< zry Q x {td > t) = for every t > 0. Thus, it suffices to show that 
lim lim E?[exp(^ TDAt )] = l. 

t|0 x—>y,x£D 

For this, note that by Jensen's inequality, 

1 > E^[exp(A TDAt )] > exp(E^[A TDAt ]). 

On the other hand, by the Cauchy-Schwarz inequality, 

(E$[A TDAt }) 2 < E x [H 2 TdM ]E x [A 2 TdM ], 

where H is the martingale given by (2.9). We know from [8] that 
sup^g^, sup sg [ 01 ] E x [ff|] < oo. Observe that 

At = N?+ C q{X s )ds = v{X t )-v(X Q )- (\v{X s )dM s + f q{X s )ds. 
JO Jo Jo 
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Since v is bounded and continuous, it is known (see, e.g., [8]) that 

-v(x). 



lim B x [(v(X TdM) 

x—ty,x£L) 



v(x) and lim B x [v(X TDAt ] 

x^y,x£D 



Thus, 



lim V x [A 2 TdM ] 

x—>y,x£D " 

< lim cE 2 



< c lim E 3 



< c lim E 3 

x^>y,x£D 



{v(X tdM 

•T D M 



v{Xo)Y 

/ rT D At 

+ y o v^(js: s )i 2 d S + ^ \ q (x s )\ 



ds 



\Vv(X s )\ 2 ds 



+ 



fr D At f rr D /\t \ 

2 J \<l( X ')\{J \q(X r )\drjds 



T D At 



\Vv(X s )\ 2 ds 



+ 2[swpG D \q\(z)\) / \q(X s )\ds 



0. 



In the second to the last inequality we used the Markov property of X, while 
in the last equality we used the fact that \Vv\ 2 and q are in the Kato class and 
the dominated convergence theorem. This proves that lim^o ^ m x^y,xeD x 
[A TDAt ] = 0. Consequently lim^ohm^^gDE^exp^^/u)] = 1. This proves 
the theorem. □ 

4. Without Markovian assumption (1.7). In this section, we will drop 
the Markovian condition (1.7) and give a general result on the existence and 
uniqueness of the weak solutions of the Dirichlet boundary value problem 
(1.8). 

Let h = (h\(x), . . . , h n (x)) :M. n — > W 1 be a measurable function such that 
h € L p (K n — > R n ) for some p> n. Let ji be a signed measure in K n . Consider 







2 . z -^ 1 dxi 

1,3=1 1 











i=l 



dxi 



The quadratic from (C, V(C)) associated with Q is given by V(C) = W^ 2 {R n ) 
and for u, v E W 1 ' 2 ^™), 



C(u,v) = (-Qu,v) 
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We will regard (C, T^ 1,2 (IR n )) as a lower-order perturbation of the symmet- 
ric Dirichlet form (£, W 1 ' 2 ^™)) associated with the infinitesimal generator 
^V(-AV). Recall that X is the symmetric diffusion process with infinitesimal 
generator ±V(AV), or equivalently, with (£, W^QR™)). Let A% be the CAF 
of X whose Revuz measure is fj,. It is proved in [8] (see also [20]) that the 
semigroup {T t ,t > 0} associated with (C, W ,2 (R n )) is given by 

rt 



T t g{x) = E 



g(X t )exp(J (A~ 1 h)(X s )dM s 



(4.2) -\ (\A- 1 \ l *{X s )ds + A<t 

2 Jo 

= Er^)e A '], 

where E^* stands for the expectation with respect to the diffusion measure 
{P* x ,x£W n }, defined by 

(4.3) 



Tt 



dP 2 

where 

(4.4) H t = exp (j\A~ 1 h)(X s ) dM s - i jT hA" 1 ^^) da) . 

Let ^ be a positive Radon measure of finite energy with respect to the sym- 
metric Dirichlet form {£ , W 1,2 (M. n )). Note that due to the L p -integrability 
of h and the Kato class condition on fx, there exists a constant ao > 1 such 
that for every a > a , 

C a (u, v) := C(u, v) + a(u, v) 

is a positive definite quadratic form satisfying 

c~ 1 C Q ,(u,u) < £\(u,u) < c a C a (u,u) for every u E W 1,2 (R n ). 

By Lax-Milgram theorem, for any a > ao, there is a unique function in 
D{£), denoted by U a u, such that 

(4.5) C a (U a u,v) = / v(x)v{dx) 



for all v G D{£). U a v is called the a-potential of v associated with the 
quadratic form (C a ,D(C)). We have the following useful representation for 

UryU. 



TIME-REVERSAL AND ELLIPTIC BOUNDARY VALUE PROBLEMS 



21 



Lemma 4.1. Let A v be the positive CAF of X associated with the smooth 
measure v in the strict sense. Then, for sufficiently large a, 



U a u(x) = E x 



Rt dAt 



E 



p 



Jo 



where 



(4.6) R { t a) =exp(^j\A~ 1 h){X s )dM s - h^ _1 h*(X s ) ds - at + A% 



Proof. By the same proof of [14], Lemma 2.2.5, we see that there exists 
an £-nest consisting of an increasing sequence {Fk,k > 1} of compact sets 
of W 1 such that \\U a (lF k v)\\oo < oo for each k > 1. By restricting to if 
necessary, we may assume that g(x) := U a u(x) is bounded in the sequel. By 
(4.5), it follows that for every v € Vr 1,2 (R n ), 



1 n f 



^ dg dv 
dxi dxj 



dx 



V / hi(x)— -v{x)dx+ [ g(x)v(x)fi(da 
jRn dxi Jr» 

+ / v{x)v(dx) — a \ g(x)v(x)dx. 



(4.7) 



By [14], Theorem 5.4.2, (4.7) implies that P z -a.s. 

rt n rt 



g(X t )=g(X )+ rVg(X s )dM s -J2 f ' h i {X s )^-{X s )ds 

JO ~{ JO OXi 

- [ g(X s )dA% + a f g(X s )ds-Al 
Jo Jo 

Note that satisfies 

dR^ = R { t a) (A^h^Xt) dM t + R { t a) dA? - aR ( t a) dt. 

By Ito's formula, it follows that 

(4.8) g(x t )Rl a) = g(X ) + K t - f* R^A", 

Jo 



where 



K t := f g(X s )R<f\A- 1 h)(X s )dM s + f R^Vg(X s )dM s 
Jo Jo 
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is a local martingale. Let {rk,k > 1} be an increasing sequence of stopping 
times approaching to oo such that {Kf ATk , t > 0} is a martingale for every 
k > 1. It follows from (4.8) that 



g{x) = E a 



t/\T k 



R{ a) dA\ 



+ B x [g(X tATn )R {a] 



thT k \- 



Note that since \n\ € K n , it follows from Khasminskii's inequality, Cauchy- 
Schwarz inequality and [8], Theorem 3.1, < s < T} is uniformly inte- 

grable under for every x € W 1 and T > 0. Letting /c — > oo, by dominated 
and monotone convergence theorems we get 



(4.9) 



g{x) = e x 



R{ a) dA\ 



+ E x [g(X t )R i 



Hi 



Since the Revuz measure [i of is in the Kato class, for sufficiently large 
a, we have 

lim E.^^O. 

t— »oo 

As (7 is assumed to be bounded, letting t — > oo in (4.9) we obtain that 



g{x) = E a 
which completes the proof. □ 



R{ a) A v Q 



Assume from now on that D is a bounded Lipschitz domain. In partic- 
ular, the boundary of D is regular with respect to (^A,D). Let R be the 
multiplicative functional of X defined as in (4.6) with a = 0; that is, 

Rt = exp^V" 1 h)(^ s ) dM s -± J h^rVpQ ds + . 

The following is a gauge theorem involving the Girsanov as well as Feynman- 
Kac functional, which extends the corresponding result in [15] where the 
diffusion matrix A is the identity matrix and the Kato class measure fi is 
absolutely continuous with respect to the Lebesgue measure [i.e., fi(dx) = 
q{x) dx]. 



Theorem 4.2. Let D be a bounded Lipschitz domain in W 1 , h e L p (R n — 
W 1 ) for some p> n and fi a signed measure in K n . Define u(x) := Ei x [R Td ]. 
Then if u(x) < oo for some x 6 D, then u is bounded between two positive 
constants in D. 
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Proof. Let Gd denote the Green function of the symmetric diffusion 
X in D and fi + , fi~ be the positive, negative part of fi, respectively. Since 
[i is in the Kato class, we have by Jensen's inequality 

inf u(x) 



> inf exp E 

~ x&D V 



[ TD (A- l h){X s ) dM s - \ r hA~ l Yi*(X s ) ds + M 
Jo 2 Jo 

A^-y^hA'^iX^ds 
> exp(-||Gz)/i"|| 00 - c||Gzj|/i| 2 ||oo)- 



inf exp E 7 

xdD * 



Let Y be the diffusion process with infinitesimal generator | J2?,j=i J5~( a «i ( x ) x 
af - ) + Ya=i ^ * s known (see [8]) that Y can be obtained from X 

through Girsanov transform by the exponential martingale H in (4.4). That 
is, if for every x 6 W 1 , we let P* be the measure defined by 



(4.10) 



dPl 



dP x 



J", 



Ht for every t > 0, 



then the diffusion process X under P* has the same distribution as Y start- 
ing from x. 

Note that R TD = H TD exp(A^ D ). For every k > 1, we have by (62) of 
Sharpe [24], 



Er>Aexp«J]=E 



Passing k 
(4.11) 



£° k A exp«) d(-H s ) + (k A exp«J)H. 

= E x [(kAexp(A? D ))H TD }. 
oo, we have 
Ef[exp(^ D )] =E x [exp(M)H 7 »] = E X [R 



~i~D 



Since h £ L P (D; dx) for some p > n, by Ancona [1], the Green function of 
Y in D is comparable to that of X in D. Hence /i is in the Kato class of Y 
in D. By [5], Theorem 2.2, we have x i— > E^* [exp(A^ D )] is either bounded 
or identically infinite on D. It then follows from (4.11) that u(x) is either 
bounded on D or identically infinite, which proves the theorem. □ 

Now consider the Dirichlet boundary value problem 

(4.12) Gu = mD with u = f ondD, 

where / € C(dD). Recall that {P^,x £ D} are the probability measures 
defined by (4.10). 
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Theorem 4.3. Recall that R is the multiplicative functional of X de- 
fined as in (4-6) with a = 0. Assume that ~E X0 [R TD ] < oo for some xq € D. 
Then there exists a unique, continuous weak solution to the Dirichlet bound- 
ary value problem (4-12), which is given by 

(4.13) u(x)=E x [R TD f(X TD )]=E**[e A ?nf(X TD )], x € D. 

Proof. We first show that u defined by (4.13) is a weak solution to the 
Dirichlet boundary value problem (4.12). Put 



v 1 (x) = E**[f(X TD )} and v 2 {x) =E 



TD 



vx{X s )e A ° dA» 



Then ui(x) = v\(x) + v 2 (x). By [8], Theorem 4.5, we know that v\ is locally 
in W 1 - 2 (M n ) and 



C(«i,0) 



1 2 

vi(x)4>{x)fj 1 (dx) for every <j) G W ' (D). 



D 



Here the quadratic form (C, W ' (D)) is defined by (4.1). We will show now 



that v 2 £ Wq' 2 {D). Let 

G (3 g(x)=Br 



TD 



e-^e A "g(X s )ds 



Denote by Gg the adjoint operator of Gg in L 2 (D;dx). Applying Lemma 
4.1 to the process X killed upon leaving D, we have 



v 2 (x) - PG L pV 2 (x) =E£ 



p 



TD 



e-^e^v^X^dA? 



Ug(vin), 



where Ug(vin) is defined as in the proof of Lemma 4.1 but with the killed 
process X D in place of X. Since /i belongs to the Kato class of X (and of 
Y by the proof of Theorem 4.2) in D, we have 



(3(v 2 ~ f3GgV 2 ,V 2 ) 

= /3 [ eH f D e-^e^v^Xt) dA* 
Jd [Jo 



v 2 (x) dx 



/3 / U 3 (vifi)v 2 (x)dx 

D 

/3G g v 2 (x)vi(x)[i(dx) 



D 



\Vvi\ 2 dx + 



D 



\v\ I 2 dx 



+ C 2 / \f3G 3 v 2 \ 2 dx + C(PGpv 2 ,(3G 3 v 2 ) 



D 



TIME-REVERSAL AND ELLIPTIC BOUNDARY VALUE PROBLEMS 25 

<Cx( [ \X7 Vl \ 2 dx+ ( \ Vl \ 2 dx) + C 2 -^ [ \l3Gnv 2 \ 2 dx 
\Jd Jd J j3-a J D 

+ \P( V 2 -PGpV2,V 2 ), 

where in the last inequality [19], Lemma 3.1 (i) , is used. Thus, 

SUP/3(W2 — f3G/3V2,V2) < oo, 
/3>0 

which implies that v 2 € W ' (D). Moreover for 4> £ Wq' 2 (D), 

C(v2,4>)= lim (3 (v2 - {3GpV2,(j>) = hm / f3Gp(p(x)v 1 (x)iJ.(dx) 

vi(x)(p(x)fi(dx). 

D 

The last equation follows from the fact that (3Gp(j){x) converges to 4> in the 
Dirichlet space W 1,2 (£>)). Thus, for <j> G Wq' 2 (D), 

C(u,4>) =C(v 1 ,</>) +C(v 2 ,<t>) =0. 

This means that u is a weak solution to Qu = in D. It is a well-known fact 
in the theory of PDE (cf. [16]) that u is continuous inside D. Next we show 
that the boundary condition is fulfilled, that is, 

(4.14) lim ui{x)=f{y) 

x—>y,x£D 

for every y £ dD. It is proved in [8] that every point in dD is a regular point 
of D c with respect to Y 

lim Vl (x) = lim Ef[/(X r J] = f(y) 

x—>y,x£L) x—>y,x£D 

for every y € dD. As u\ = v\ + V2, it suffices to show that 

lim V2(x) = for y £ dD. 
x^y,x&D 

Note that 

V2{x)=vT[f{X TD ){e A *D-l)]. 
For any t > 0, write V2 as 

U2(x) = Ef [f(X TD )(e A "o - l)- TD <t} + Ef [/(X r J( A - 1) ;TD > t]. 
Now 

[Ef[/(X r J( e ^ -1);td>*]| < ll/IUEfKA +l);r I? >t] 

= ||/|| 00 Er[e^EP;[e^+l];r D >t] 

<||/|| 00 Er[e^ ! (M + l);r D >t]. 
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Here in the last inequality, we used the assumption that E x [i? TD ] < oo 
for some x € D and thus by Theorem 4.2, M := sup^^E^ [exp(j4^ D ] = 
sup^g^, ~E X [R TD \ < oo. Let y E dD. Since 

lim P* (t d > t) = and sup E^* [e A ^ ] < oo, 

x^y,x&D xG £) 



we have 



lim EP [/(X T J(e A ^-l):r D >t]=0. 

x— >y,x£V 



Thus for every t > 0, 

lim |«2(x)| = limsup \E^[f{X TD ){e A "n - \)- TD < t}\ 

x^y,xeD x^y,xeD 

< H/lloo limsup E**[e A ^ -l;r D <t] 

x—>y,x£D 

< H/IU limsup (E x [(F i ) 2 ]) 1 / 2 (E x [(e A i Ml - 1) 2 ]) 1/2 . 

x^y,x£D 

By [8], Theorem 3.1, we know that sup xg£ ) E x [(H t ) 2 ] < oo, where H is de- 
fined in (4.10). On the other hand, since \x is in the Kato class of X, it 
follows from the Khasminskii's inequality that 

lim sup E x [(e A *"' - l) 2 ] = 0. 

Thus we conclude that lim x ^ ytXI =D V2(x) = 0. So (4.14) is established. 

To prove the uniqueness, let u\ be any, bounded continuous weak solution 
of the Dirichlet boundary value problem (4.12). Then for <f) G Wq ,2 (D), 

du\{x) 



^hiix) — 4>(x)dx 



(4.15) + \_u\{x)<l)(x)n{dx) 

du\{x 



f au\[x) f 

V / hi(x)— 4>{x)dx+ / ui(x)(j)(x)n{dx) 

~[Jd oxi Jd 



By [14], Theorem 5.4.2, it follows from (4.15) that the following decomposi- 
tion holds: 



(X t )- Ul (X ) = / V Ul (X s )dM s -y2 / hi(X s )-±(X s )ds 

JO ^ Jo OXi 

Ul (X s )dA% 
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for t <td- Recall that 

Rt = exp^V" 1 h)(^ s ) dM s - 1 jT hA-^X,) ds + A? 
satisfies 

dR t = i?i(^ _1 h)(X t ) dM t + R t dA^. 

Applying Ito's formula, we get that for t<Tjj, 

(4.16) d(u 1 (Xt)Rt) = RtVu 1 (Xt)dMt + u 1 (X t )Rt(A- 1 h)(Xt)dM t . 

This shows that {u\{Xt^ TD )Rt/\ TD ,t > 0} is a P^-local martingale for every 
x € D. We claim that {R t /\T D ,t > 0} is uniformly integrable with respect to 
P x for every x € D. To see this, write 

Rt/\T D = RT D l{t>T D } + Rtl{t<T D }- 

Obviously the first term {R TD l^ t > TD }} is uniformly integrable. For the sec- 
ond term, by Jensen's inequality, 

1 {t<T D }^x[Rr D \^t\ 

= l {t<TD} E x [Rt(R TD o9 t )\T t ] 

= Rti{t<T D }Ex t [Rt d ] 





sup E x 






\l 


x£D 


Jo 



hA- 1 h*(X a )ds 



>cRt±{t <TD } 

for some positive constant c, where we have used the fact that \h\ 2 , \i are both 
in the Kato class. The above inequality implies that {Rtl{t<r D }-> t > 0} is 
also P^-uniformly integrable. Therefore, {RtAr D , t > 0} is P x -uniformly inte- 
grable. Since u is bounded continuous, we conclude that {u\ {Xt/\ TD )Rthr D , t > 
0} is uniformly integrable, hence a P x -martingale for every x £ D. Conse- 
quently, 

E x [ui(X tA T D )RtAT D ] =ui(x). 
Letting t — > oo, we get that 

Ul (x)=E x [f(X TD )R TD ], 
which proves the uniqueness. □ 

Now we can get to the main results of this section. Define 

Z t = exp j^V -1 b)(A%) dM s + ^(A^bXXs) dM^j o r t 



(4.17) 



X - J\h - h)A~\h - b)*(X s ) ds + J\{X S ) ds\. 
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Recall that X is the symmetric diffusion with infinitesimal generator 
\V(AV) and M is the martingale part of X in (1.9). For domain D C W 1 , 
tjj := inf{i > 0:Xt D} is the first exit time from D by diffusion X. The 
following theorem is a new type of gauge theorem, in comparison with those 
found in [2, 5, 9]. 

Theorem 4.4. Let D be a bounded Lipschitz domain contained in some 
ball B R , A be an n x n symmetric positive definitive matrix satisfying the 
condition (3.1) of Lemma 3.1, |b| + |b| G L p (D;dx) for some p> n, and 
IdQ G K n . Then Z TD of (4-17) is well defined under P x for every x G D. If 
~E XQ [Z TD ] < oo for some xq G D, then the function x \—> ~E X [Z TD ] is bounded 
between two positive constants on D. 

Proof. As before, put 

M t = /*(A~ 1 b)(X s ) dM s for t > 0. 

Let R > so that D C Br := B(Q, R). By Lemma 3.2, there exits a bounded 
function v G wl' p {B R ) C Wq' 2 {B r ) such that 

M t or t = -M t + 

and that v G Wq 1,2 {Br) satisfies the following equation in the distributional 
sense: 

(4.18) div(AVv) = -2div(b) in B R . 

Note that by Sobolev embedding theorem, v G C(M n ) if we extend v = 
on D c . Thus both M and iV" are CAFs of X in the strict sense (cf. [13], 
Theorem 1), and so is 1 1— > Mt o rt- Moreover, 

^ D (A~ i: b)(X s )dM s ^ or TD 

= ~ (A- 1 h){X s )dM s + N V TD 

= - r {A~ 1 h)(X s )dM s + v{X TD ) 
Jo 

= - r (A~ i: b) (X s ) dM s + v{X TD ) 
Jo 

- ( TD Vv(X s )dM s . 
Jo 



- v(X ) - M\ 

- v(X ) 
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Thus 

V(X TD ) / ,-T D 

= ~Kx^ e M I ( A ~ ( b - b ) - Vv )( X s) dM s 

(4.19) ° 

+ jT^-^b-bM-^b-br)^)^). 

So Z TD is well denned under for every x G D. Since t> is bounded, b — 
b G L p (D;dx), Vu € L p (D;dx) and l^g € K n , the theorem follows from 
Theorem 4.2. □ 

Recall that C is the second-order differential operator defined by (1.3). 

Theorem 4.5. Let D be a bounded Lipschitz domain contained in some 
ball Br, A be an n x n symmetric positive definitive matrix satisfying the 
condition (3.1) of Lemma 3.1, |b| + |b| G L p (D;dx) for some p> n, and 
IdQ £ K n . Let Z be defined in (4-17) and assume that E X [Z TD ] < oo for 
some x G D. Then for every f G C(dD), there exists a unique weak solution 
u to Cu = in D that is continuous on D with u = f on dD. Moreover, the 
solution u admits the following representation: 

(4.20) u(x)=E x [Z TD f(X TD )] forxeD. 

Proof. Let u be defined by the right-hand side of (4.20). Recall that v 
is the function in Wq ,p {Br) that is continuous on W 1 if we extend v = off 
D in the proof of Theorem 4.4. Define 

Z TD := exp^^ (A-^b - b - AVv))(X s ) dM s 

--/ (b-b- AVv)A" 1 (h-b- AVv)*(X s )ds 
2 Jo 

- [ TD {b-b,Vv)(X s )ds 
Jo 

+ -y o {Vv)A(Vv)*(X s )ds + j o q(X s )ds). 

We have by (4.19) that 

Ztd = e v(X a ) Zt o- 

Since the function v is bounded and continuous on D, it follows that E XQ [Z TD ] < 
+oo for some xq G D if and only if E Xo [Z TD ] < +oo. Define g := e v u. Then 

(4.21) g(x)=E x [Z TD (e v f)(X TD )} for x G D. 



30 Z.-Q. CHEN AND T. ZHANG 

Let 

1,3=1 1 x J ' 

n o 

- (b - b, Vv)(x) + ^(Vv)A(X7v)*(x)+q(x). 

By Theorem 4.3, g is a weak solution to the Dirichlet boundary value prob- 
lem 

C\g = in D and g = fe v on dD. 

Moreover, g is continuous on D with g = fe v on dD. Hence, u = e~ v g is 
bounded and continuous on D with u = f on dD. Note that for any ip € 
Wt> 2 (D), 

Q*{g^):={-C l9 ^) 

_if /• , 

~ 2 i A±J D aij dx i dx j dX 



(4.22) -J2 f ^ ^-(AVv^^dx 

~[ JD OXi 

— / qgijjdx+ / (b — b, Vv)g%j) dx 

JD il" 



1 

~ 2 
0. 



Recalling the definition of the quadratic form Q from (1.6), we thus have 
for every <j> £ C™(D), 

2 r-^jD dxi dxj 'HJd dx { 



*>j'=i' 



cfidx 



(4.23) 



— / bi— — ge v dx — [ qge v 
Jd oxi Jd 

1 " f dg 9(e~» 1 " /■ dgde-" 

2 j D ^dx--dx— dX ~ 2 .ln"' J d7,^-"' Lr 
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1 " f 8e~ v 86 n f dg 

2, J — 1 % — 1 

n f 8e~ v n f ~ 86 

b ^ — 94>dx-Y] / bi—ge~ v dx 
~{JD OXi i^i-' D dx i 

qge~ v (j)dx. 



D 



Applying (4.22) with tp = 6e " we obtain 

if / a ( x ) gg d ^~ V ) dx 
2^L a ^ X) d Xi dx, dX 



= J2 I (bi(x) -bi(x) - (aVvUx))^-6e- v dx 

jr{ Jw- axi 

+ / q(x)g(x)6e~ v dx 

- f (b - h,Vv){x)g{x)6e- v dx 

+ -/ (Vv)A(Vv)*(x)g6e- v dx. 

2 JR" 

Substituting this expression into (4.23), we get after cancellations that 

MgS) 



Q(u,6) = -J2[ h(x)^e^dx-J2 f (AVv)i( 
Jm" oxi ~[ m n 



x)— — 6e v dx 

OXi 



+ / {h,Vv)6ge~ v dx + - f (Vv)A(Vv)* g6e~ v dx 

(4.24) 

1 n f dg 8e~ v 



1 A f , ,de~ v 86 , 
2,, ^R" <9xj dec,- 



»,j'=i ' 



In the sequel, we write div(-) for the divergence in the distribution sense. 
Now, 

OXi 



(4.25) = / div(be~ v )g6dx 
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= / div(b)e~ v g</>dx- [ (b, Vv)<pge~ v dx. 
By virtue of (4.18), 

/ dw{b)e- v g(f>dx 

= -- ( dw(AVv)e~ v g(f)dx 

(4.26) =- f (AVv,V{e- v g<p))dx-- f {AVv,Vv)e~ v gcj)dx 

2 Jw 1 2 it" 

+ \f {AVv,Wg)e- v ct>)dx 
+ - / (AVv,V(/))e- v g4>dx. 

2 jRn 

Combining (4.24)-(4.26), we arrive at 

Q(u, </>) = -/ (A(x)Wv(x) ■ Vg(x))0(x)e- v{x) dx 

JR™ 

1 n f da de~ v 1 f 

"oE / ^ji^j^^^^^+n (AVv,Vg)e- v (l)dx 

= 0. 

This shows that u is a weak solution to Cu = in D. Recall that we have 
showed earlier that u is continuous on D with u = f on 3D. 

It remains to show the uniqueness. Suppose that u is a weak solution to 
Cu = in D and that u is bounded and continuous on D with u = f on 
dD. Then, for any G W 1,2 (D), Q(u,$ = 0. Let u be the solution of (4.18). 
Put g(x) = e v ^u(x). Running the above proof backward, we can show that 

1 2 

Q*{g,ip) = for any V € W c \' (D). Hence g is a weak solution to Cig = 
in D that is bounded and continuous on D with g = fe v on 3D. It follows 
from Theorem 4.3 that g is given by (4.21). Arguing similarly as at the 
beginning of the proof, we conclude that u can be expressed as in (4.20). 
This establishes the uniqueness. □ 

Remark 4.6. For b 6 L p (M n ; dx), by [26], Theorem 4.1, there is a unique 
bounded weak solution v € Wq' 2 (Bji) such that 

div(^Vt;) = -2div(b) in B R . 

Condition (3.1) of Lemma 3.1 is used to guarantee that Vf € L p {Bn;dx) 
for some p > n [see (4.18)], which in turn by the Sobolev embedding the- 
orem implies that v € C(R n ) if we extend v = on D c . Condition (3.1) 



TIME-REVERSAL AND ELLIPTIC BOUNDARY VALUE PROBLEMS 33 

can be dropped from Theorems 4.4 and 4.5 if we assume directly that 
Vi> € L p {Br; dx) for some p> n. 

5. Generalized Feynman Kac transform. In this section, we consider the 
special case of (1.3) with b = b = — AVp, for some p € W 1,2 (M. n ) with |Vp| 2 £ 
L p (R n ;dx) for some p> n. Note that by Sobolev embedding theorem, p € 
C(K n ). As mentioned in the Introduction, the quadratic form (Q, W 1 ' 2 (W n )) 
in (1.6) takes the following form: 

. 1 A f . . du dv , 1 A / , .d(uv) dp 
2 ^ 7 R ~ dxi dxj 2 f-±. jRn dxi dxj 

(5.1) 

+ / u(x)v{x)q(x) dx 



dx 



for u,v € W 1 ' 2 (W l ). Let X be the symmetric diffusion with infinitesimal 
generator ^V(AV). When q = 0, it is established in Chen and Zhang [7] 
that (Q, W 1 ' 2 (W 1 )) is the quadratic form associated with the generalized 
Feynman-Kac semigroup {T t ,t > 0} defined by 

T t f(x)=E x [e N "f(X t )}, x£R n ,t>0, 

where N p is the CAF of X in the Fukushima's decomposition of 

p(X t )-p(X ) = M t p + N t p , t>0. 

(When A is the identity matrix, the above result is proved in [17].) Recall 
that M is the martingale defined in (1.9) and that Mf = J v p(X s ) dM s . 
Note that (cf. [3]) 

N t p = -±(M t p + M t p or t ), t>0. 

We thus have from (4.17) that, under the condition of Theorem 4.5, for every 
/ G C(8D), 



u(x) :=B X 



exp^ D + q(X s )ds^f(X TD : 



xedD, 



is the unique weak solution to Cu = in D that is continuous on D with 
u = f on dD. However, we can do better in this case. More specifically, we 
can drop condition (3.1) of Lemma 3.1 in this case. 

Theorem 5.1. Let D be a bounded Lipschitz domain in W l . Let A be 
a symmetric n x n bounded positive definite measurable matrix on W 1 , p 6 
W ,2 (K. n ) with Vp G L p (W n ;dx) for some p> n. Let operator C be defined 
by (1.3) with b = b = —AVp and q £ K n . Define 

Z t = exp(N p + j\(X s )d s y t>0. 
Then Z-t is well defined under P x for every x G W 1 . 
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(i) If ~E X0 [Z TD ] < oo for some xq E D, then x i— ► E X [Z TD ] is bounded 
between two positive constants on D. 

(ii) Suppose that E XQ [Z TD ] < oo for some xo £ D and f E C(dD). Then 

u(x):=E x [Z TD f(X TD )}, xED, 

is the unique weak solution of Cu = in D that is continuous on D with 
u = f on dD . 

Proof, (i) Since Vp £ L p (IR™;<ix) for some p> n, by Sobolev embed- 
ding theorem, p can be taken to be continuous on W 1 . By [13], Theorem 1, 
the following Fukushima decomposition in the strict sense holds: 

p(X t )-p(X ) = M t p + N t p , t>0, 

where M p is the MAF of X in the strict sense of finite energy and N p is the 
CAF of X in the strict sense of zero energy. It follows 

Z TD = exp (p(X TD ) - p(X ) - MP D + J T J q(X s ) ds 

exp(p(X r J) 
exp(n(X )) TD 

is well defined under P x for every x E W 1 , where 

R t :=exJ- [\p{X s )dM s - \ f\(Vp)A(Vp)*)(X s )ds 
(5.2) W ° 2J ° 

+ J*fyvp)A(Vp)* + q yx s )ds 

Since p is bounded on D and Vp £ L p (E. n ; dx), E :ro [Z T£) ] < oo if and only if 
Ea; [-R TD ] < oo, and so the conclusion of (i) follows from Theorem 4.2. 
(ii) By (5.2), 

e p ^u(x) = B x [R TD (e p f)(X TD )} for x E D. 

Define 

^E^r KW)^: - A(x)(Vp(x)y ■ V + - (Vp)A(Vp)* + ?■ 



It follows from Theorem 4.3 that v := e p u is the unique weak solution for 
Qv = in D that is continuous on D with v = e p f on Unwinding it for 
u similar to that in the proof of Theorem 4.5, we conclude that u £ W^' (D) 
is the unique weak solution of Cu = in D, that continuous on D with u = f 
on dD. □ 
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Remark 5.2. To better illustrate the main ideas of this paper, we have 
not strived to establish our theorems in the most general possible form. 
For example, in Theorems 4.4, 4.5 and 5.1, the potential function q can be 
replaced by a signed Kato class measure /i, just like in Theorems 4.2 and 
4.3. 

Acknowledgments. We thank the anonymous referee and K. Kuwae for 
helpful comments on an earlier version of this paper. 
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